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Abstract. The primary goal of the GERmanium Detector Array (Gerda) 
experiment at the Laboratori Nazionali del Gran Sasso of INFN is the search 
for the neutrinolcss double beta decay of 76 Ge. High-purity germanium detectors 
made from material enriched in 76 Ge are operated directly immersed in liquid 
argon, allowing for a substantial reduction of the background with respect to 
predecessor experiments. The first 5.04 kg-yr of data collected in Phase I of 
the experiment have been analyzed to measure the half-life of the neutrino- 
accompanied double beta decay of 76 Ge. The observed spectrum in the energy 
range between 600 and 1800 keV is dominated by the double beta decay of 76 Ge. 
The half-life extracted from Gerda data is if" = (1.841^;^) • 10 21 yr. 
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1. Introduction and scope 

Neutrinoless double beta (0^/3/3) decay of atomic nuclei (A,Z) —> (A.Z + 2) +2e~ is a 
forbidden process in the Standard Model (SM) of particle physics because it violates 
lepton number by two units. An observation of such a decay would demonstrate 
lepton number violation in nature and would prove that neutrinos have a Majorana 
component. For recent reviews, see [1]. The two-neutrino double beta (2^/3/3) decay 
of atomic nuclei, 

(A, Z) (A, Z + 2) + 2e~ + 2v e , 

with the simultaneous emission of two electrons and two anti-neutrinos, conserves lep- 
ton number and is allowed within the SM, independent of the nature of the neutrino. 
Being a higher-order process, it is characterized by an extremely low decay rate: so 
far it is the rarest decay observed in laboratory experiments. It is observable for a 
few even-even nuclei and was detected to-date for eleven nuclides; the corresponding 
half-lives are in the range of 7 ■ 10 18 - 2 ■ 10 24 yr [2, 3, 4]. 

The measurement of the half-life of the 2^/3/3 decay (T 2 ^) is of substantial in- 
terest. For example, model predictions of the 0^/3/3 half-life require the evaluation 
of nuclear matrix elements. These calculations are complicated and have large un- 
certainties. They are different from those required for the 2v[3(3 decay, but it has 
been suggested [5, 6] that, within the same model framework, some constraints on the 
Of/3/3 matrix elements NME°" can be derived from the knowledge of the 2v(5j3 nuclear 
matrix elements NME 2 ". Also, the nuclear matrix element NME 2 " which is extracted 
from the measurement of the half-life of the 2^/3/3 decay can be directly compared 
with the predictions based on charge exchange experiments [7, 8]. A good agreement 
would indicate that the reaction mechanisms and the nuclear structure aspects that 
are involved in the 2^/3/3 decay are well understood. 

The GERmanium Detector Array (Gerda) experiment at the Laboratori Nazion- 
ali del Gran Sasso of INFN searches for the 0^/3/3 decay of 76 Ge. High-purity germa- 
nium detectors isotopically enriched in 76 Ge are operated bare and immersed in liquid 
argon in order to greatly reduce the environmental background. As a first result of 
this ongoing research, the present paper reports a precise measurement of the half-life 
of the 2z//3/3 decay of 76 Ge. The data used in this work encompasses an exposure of 
5.04 kg-yr, taken between November 2011 and March 2012. 



2. The Gerda experimental setup 

A brief outline of the components of the Gerda detector that are most relevant for 
this work is given below; a detailed description can be found in [9, 10]. 

The Gerda experimental setup is shown in figure 1. At the core of the setup 
there is an array of high-purity germanium detectors (HPGe) . They are operated bare 
in liquid argon (LAr) which acts both as a coolant and as a shield against the resid- 
ual environmental background. The array configuration consists of eleven germanium 
detectors: eight are made from isotopically modified germanium ( cnr Ge), enriched to 
about 86% in 76 Ge, and three are made from natural germanium ( nat Ge), with a total 
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Figure 1. Artist's view of the Gerda experiment. The detector array is not to 
scale. 



mass of 17.67 kg and 7.59 kg, respectively. The enriched detectors come from the 
former Heidelberg-Moscow (HdM) [11] and Igex [12] experiments. They underwent 
specific refurbishing processes before operation in Gerda [13, 14]. The germanium 
detectors are mounted in strings with typically three diodes each. Signals are ampli- 
fied by low noise, low radioactivity charge sensitive preamplifiers [15] with 30 MHz 
bandwidth operated inside the LAr. They are digitized by a 14-bit 100 MHz contin- 
uously running ADC (FADC) equipped with anti-aliasing bandwidth filters. In the 
offline analysis the waveforms are digitally processed to reconstruct the event energy. 

The detector array is surrounded by 64 m 3 of 5.0-grade LAr, contained in a vac- 
uum insulated cryostat made of stainless steel, lined on the inner side by a 3 to 6 cm 
thick layer of copper. The cryostat is in turn placed at the centre of a 580 m 3 volume 
of ultra-pure water equipped with 66 photomultiplier tubes to veto the residual cosmic 
ray muons by the detection of Cherenkov light. The large water volume also serves as 
a shield to moderate and capture neutrons produced by natural radioactivity and in 
muon-induced hadronic showers. 

The energy scale is set by using calibration curves, parametrized as second-order 
polynomials, derived for each detector by calibration runs taken with 228 Th sources. 
The stability of the energy scale is monitored by performing such calibration runs ev- 
ery one or two weeks. Moreover, the stability of the system is continuously monitored 
by using ad hoc charge pulses generated by a spectroscopy pulser that are regularly 
injected in the input of the charge sensitive preamplifier. 

All Gerda detectors but two exhibit a reverse current of the order of tens of 
pA. The two problematic detectors showed an increase of leakage current soon after 
the beginning of their operation; therefore their bias high voltage had to be reduced 
and finally completely removed. These two detectors do not contribute to the present 
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data set, although they are accounted for in the detector anti-coincidence cut de- 
scribed below. The total mass of the operational enriched detectors is 14.63 kg. The 
nat Ge detectors are not considered in the present analysis because of their low content 
of 76 Ge. The average energy resolution for the six enriched detectors considered in 
this study is about 4.5 keV (full width at half maximum) at the Q^-value of 76 Ge 
(2039 keV) [10]. 

A key parameter which is required for the computation of Ty 2 is the total 76 Ge 
active mass. It is calculated as the product of the total detector mass with the isotopic 
abundance of 76 Ge (/re) and the fraction of the active mass (fact)- 
The average 76 Ge isotopic abundance of the six enriched detectors considered in this 
work is 86.3% [10, 12, 16, 17]. 

In Gerda p-type semi-coaxial detectors are used, for which a part of the volume 
close to the outer surface is inactive. After mechanical modifications and processing 
of the germanium diodes at the manufacturer, their active mass was re-determined 
experimentally [10, 13]. The measured average fraction f act is 86.7%, with individual 
detector f act uncertainties of about 6.5%. 

3. The data set 

The data set considered for the analysis was taken between November 9, 2011, and 
March 21, 2012, for a total of 125.9 live days, amounting to an exposure of 5.04 kg-yr. 
Digitized charge pulses from the detectors are analyzed with the software tool Gela- 
tio [18] according to the Gerda standard procedure [19]. The pulse amplitude is 
reconstructed offline by applying an approximated Gaussian filter with an integration 
time of 5 fis (which is sufficient to avoid losses due to ballistic effects). Events gener- 
ated by discharges or due to electromagnetic noise are rejected by using a set of quality 
cuts. Due to the low counting rate, the data set has a negligible contamination of pile- 
up events. The combined efficiency of the trigger and of the offline data processing is 
practically 100% above 100 keV. Similarly, no loss of physical events is expected from 
the application of the quality cuts, as deduced by dedicated Monte Carlo studies and 
by the analysis of pulser data. Events that are in coincidence with a valid veto signal 
from the muon detector and events having a signal in more than one HPGe detector 
are excluded from the analysis. The time window for the coincidence between the 
muon detector and the HPGe detectors was set to 8 fj,s while that between different 
HPGe detectors was set to a few /is. Due to the very low muon flux in the Gran Sasso 
laboratory, the dead time induced by the muon veto is negligible. 

Given the half-life of the 2v(3(3 decay reported in the literature (about 1.5 -10 21 yr), 
the anticipated count rate of the enr Ge detectors is about 100 counts/day in the entire 
energy range up to Q^=2039 keV. Since the detectors are submerged in LAr, the 
radioactive decay of 39 Ar, which is a long-lived /3 emitter produced by cosmogenic 
activation of natural argon in the atmosphere, gives a large contribution up to its Qp- 
value of 565 keV. In fact, the low-energy spectrum is dominated by these /3 particles 
and their Bremsstrahlung photons, which account for about 1000 counts/day above 
100 keV. The 2v(3f3 decay is expected to be the major contributor in the energy 
spectrum above the end-point of the 39 Ar spectrum. For this reason, the analysis 
of the 2vj3f3 decay is performed in the range between 600 keV, which is comfortably 
above the end-point of the 39 Ar spectrum, and 1800 keV. The sum spectrum of the 
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Figure 2. Energy spectrum of the six enr Ge detectors. The 39 Ar /3 decay 
dominates the counting rate at low energies up to its Q^-value of 565 keV. The 
energy range which is considered for the 2v/3/3 analysis and a few prominent 7 
lines are also shown. 



six enr Ge detectors considered in this work is displayed in figure 2. The analysis range 
contains 8796 events in total. The probability for a 2v(3(i decay taking place in the 
active volume of the enr Ge detectors to produce a total energy release between 600 and 
1800 keV is about 63.5%; the energy range above 1800 keV is practically insensitive 
to the 2v[3f3 signal, as the probability for a decay to produce an energy release in this 
region is < 0.02%. These estimates are based on the Monte Carlo simulation of 2v(3(3 
decays in the Gerda detectors as described in section 4. Hence, the energy region 
chosen for the analysis is well suited for the study of the 2z//3/3 decay signal. 

4. Data analysis 

4-1. Statistical treatment and fit model 

The experimental spectra of the six diodes are analyzed following the binned maxi- 
mum likelihood approach described in [20]. The analysis region is divided into 40 bins, 
each 30 keV wide. A global model is fitted to the observed energy spectra. The model 
contains the 2vfjj3 decay of 76 Ge and three independent background contributions, 
namely 42 K, 214 Bi and 40 K. The presence of these background sources is established 
by the observation of their characteristic 7 lines: 1525 keV from 42 K; 1460 keV from 
40 K; 609 keV and 1764 keV from 214 Bi; 352 keV from 214 Pb (progenitor of 214 Bi). 

42 K is a short-lived (3 emitter (Qp= 3525 keV, T 1 / 2 = 12.6 h) which is present in 
Gerda as a progeny of the long-lived 42 Ar (Qp= 599 keV, T 1/2 = 32.9 yr). 42 Ar is a 
trace radioactive contaminant expected in natural argon and is produced by cosmo- 
genic activation. 214 Bi from the 238 U decay series and 40 K are 7 emitters from the 
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environmental radioactivity. 

A few more candidate 7 lines have been identified in the Gerda spectrum [10], which 
indicate the presence of small additional background contributions: in particular 208 Tl 
and 228 Ac from the 232 Th decay series, and 60 Co. However, most of the candidate 7 
lines have either a poor statistical significance and/or are seen in some detectors only. 
Given the lack of discriminating power in the data, the background contributions other 
than 42 K, 214 Bi and 40 K are not included in the fit. However, their possible impact 
on the extracted half-life Tjy 2 is included in the systematic uncertainty, as discussed 
in section 4.2; their cumulative contribution to the background is estimated to be of 
a few percent. 

The half-life of the 2v(3f3 decay is common in the fit to the six spectra of the 
cnr Ge detectors. The intensities of the background components are independent for 
each detector. The active mass and the 76 Ge abundance of each detector are also left 
free in the fit; they are treated as nuisance parameters and integrated over at the end 
of the analysis. 

The shapes of the energy spectra for the model components (signal and three 
backgrounds) are derived by a Monte Carlo simulation for each detector individually. 
The simulation is performed using the MaGe framework [21] based on Geant4 [22]. 
Assumptions have to be made in the Monte Carlo simulation about the location and 
the primary spectrum of each component of the model. The spectrum of the two 
electrons emitted in the 2vf3fS decay of 76 Ge is sampled according to the distribution 
of [3] that is implemented in the code decayO [23]. Electrons are propagated in the 
Gerda simulated setup by MaGe and the total energy released in the active mass of 
the enriched detectors is registered. 

The 42 K activity is uniformly distributed in the liquid argon volume. The decay prod- 
ucts of the j3 decay of 42 K are taken into account as the initial state in the simulation. 
However the energy deposit in the detectors is mainly given by the 1525 keV 7-ray 
(full energy peak and Compton continuum): the contribution due to the j3 particles 
is small (less than a few percent) with respect to the 1525 keV 7-ray. The actual 
position of the 40 K and 214 Bi emitters contributing to the Gerda spectrum is not 
known in detail: the assumption is made in the Monte Carlo of "close sources" . The 
ratio of the intensities of the 21 Bi 7 lines observed in the experimental spectrum is 
consistent with such an assumption. The impact on T!f/2 due to the lack of knowledge 
about the source position - which affects the peak-to-Compton ratio - is accounted as 
a systematic uncertainty. The effect of muon-induced events that are not accompanied 
by a veto signal, e.g. due to inefficiency, is estimated to be < 0.02% in the energy 
range of this analysis. 

The spectral fit is performed using the Bayesian Analysis Toolkit Bat [24] . A flat 
distribution between and 10 22 yr is taken as the prior probability density function 
(pdf) for Ty 2 . The prior pdf for the active mass fraction of each detector is mod- 
elled as a Gaussian distribution, having mean value and standard deviation according 
to the measurements performed in [13]. The analysis accounts for the fact that the 
uncertainties on the active masses f ac t are partially correlated, because of the experi- 
mental procedure employed for the measurement. Uncertainties on f act are split into 
an uncorrelated term — which is specific of each detector — and a common correlated 
term. The prior pdf for the 76 Ge isotopic abundance fa for each detector is also a 
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Table 1. Summary table of the six enriched detectors used in this work. For each 
detector, the total mass, the active mass and the isotopic abundance of 76 Ge are 
given. The uncertainties reported for the active masses are the uncorrelated (first) 
and correlated (second) errors. The last column reports the T^J 2 which is obtained 
by performing the analysis on each detector individually. The uncertainty on T^J 2 
is the fit uncertainty only. 



Gaussian, having mean value according to the earlier measurements reported in [10]. 
The uncertainty of /76 is between 1% and 3% for the single detectors. It was esti- 
mated from the dispersion of independent measurements performed with isotopically 
enriched material; this dispersion is much larger than the quoted uncertainty of each 
individual measurement. 

The fit has 32 free parameters: one common value for the 2v[3{3 half-life and 31 nui- 
sance parameters. The nuisance parameters are the common term which describes the 
correlated uncertainty of the active masses plus five independent variables for each of 
the six detectors (active mass, 76 Ge abundance and three background components). 
The detector parameters used for the prior pdf's are summarized in table 1. 

Figure 3 shows experimental data together with the best fit model for the sum 
of the six detectors. The analysis energy window contains 8796 events. The best fit 
model has an expectation of 8797.0 events, divided as follows: 7030.1 (79.9%) from 
the 2v/3P decay of 76 Ge; 1244.6 (14.1%) from 42 K; 335.5 (3.8%) from 214 Bi; and 186.8 
(2.1%) from 40 K. The individual components derived from the fit are also shown in 
figure 3. The signal-to-background ratio in the region 600—1800 keV is on average 
4:1, which is much better than for any past experiment which observed the 2v(3f3 
decay of 76 Ge. The best ratio achieved so far for 76 Ge was approximately 1 : 1 by 
HdM [17]. The model is able to reproduce well the experimental data, as shown in 
the lower panel of figure 3. The p- value of the fit derived from the procedure of [25], 
is p — 0.77. All nuisance parameters are eventually integrated over in order to derive 
the posterior pdf p{T^ 2 ) for the 2vf3f3 half-life. The p(T^J 2 ) distribution is nearly 
Gaussian; the best estimate of the half-life is 

r 2 ; 2 = (1.841°;°!) • 10 21 yr, (1) 

(fit error only). This uncertainty is calculated as the smallest interval containing 68% 
probability of p(T 2 J 2 ). It includes the uncertainty induced on T 2 ^ 2 by the nuisance 
parameters of the fit and accounts for parameter correlations. Active masses and 76 Ge 
isotopic abundances drive the fit uncertainty on Ty 2 : if these parameters were known 
without uncertainties, the la uncertainty from the fit would be about 0.03 • 10 21 yr. 
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Figure 3. Upper and middle panels: Experimental data (markers) and the best 
fit model (black histogram) for the sum of the six detectors together (linear and 
logarithmic scale). Individual contributions from 2vj3/3 decay (red), 42 K (blue), 
40 K (purple) and 214 Bi (green) are shown separately. The shaded band covers 
the 68% probability range for the data calculated from the expected event counts 
of the best fit model. Lower panel: ratio between experimental data and the 
prediction of the best fit model. The green, yellow and red regions are the smallest 
intervals containing 68%, 95% and 99.9% probability for the ratio assuming the 
best fit parameters, respectively [26]. 
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Table 2. Summary table of the systematic uncertainties on T^J 2 which are taken 
into account for this work and which are not included in the fitting procedure. 

As a cross-check, the same procedure is run for each detector separately. The 
resulting T 2 ^ values are summarized in table 1; they are mutually consistent within 
their uncertainties (\ 2 jv = 3.02/5). 

4-2. Systematic uncertainties 

The items which are taken into account as possible systematic uncertainties of T 2 ^ 
and which are not included in the fitting procedure are summarized in table 2. They 
can be divided into three main categories: (1) uncertainties related to the fit model 
(choice of the components, shape of input spectra); (2) uncertainties due to the Monte 
Carlo simulation regarding the precision of the geometry model and the accuracy of 
the tracking of particles; (3) uncertainties due to data acquisition and handling. The 
latter term turns out to be negligible with respect to the others. The most relevant 
items from table 2 are briefly discussed in the following. 

Additional background components that are not accounted for in the fit model 
might be present in the Gerda spectrum (sec ref. [10] for a list of the +-ray lines de- 
tected in the Gerda spectrum and of the corresponding intensities) . Due to the large 
signal-to-background ratio and the limited exposure these background components 
cannot be identified unambiguously. The uncertainty arising from such possible con- 
tributions is estimated to be +5.3%. Since any further background component would 
lead to a longer T^J 2 , this uncertainty is asymmetric. It is estimated by performing 
a fit with the contributions from 60 Co, 228 Ac, and a flat background added to the 
model. These additional components are treated in the same way as the "standard" 
background components ( 42 K, 40 K, and 214 Bi). The spectra from 60 Co and 228 Ac are 
simulated by Monte Carlo assuming close sources and one additional parameter for 
each detector and each additional background contribution is included in the fit. Also 
for the flat background an individual contribution is considered for each detector. The 
flat component describes the contribution coming from 208 T1 decays from the 232 Th 
chain: given the small number of events expected in the analysis energy window, this 
contribution can be roughly approximated to be constant. To a first approximation, 
also other possible non-identified background components can be accounted by the 
constant contribution to the model. 
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The systematic uncertainty on T^Y 2 due to the uncertainties in the spectra of the 
standard background components ( 42 K, 40 K, and 214 Bi) is estimated to be 2.1%. It 
is evaluated by repeating the analysis with different assumptions on the position and 
distribution of the sources and with artificial variations (e.g. via a scaling factor) of 
the ratio between the full-energy peaks and the Compton continua. 

The primary spectrum of the 2vf3/3 decay which is fed into the Monte Carlo simu- 
lation is generated by the code decayO. DecayO implements the algorithm described 
in [3], which is based on [27, 28]. The 2^/3/3 decay distributions of [3] are in principle 
more precise than those based on the Primakoff- Rosen approximation [29] . They have 
been cross-checked against the high-statistics data of the Nemo experiment for several 
nuclei: 82 Se, 96 Zn and 150 Nd [30]. The Ivfifi spectrum derived by the Primakoff- Rosen 
approximation was used in earlier works with 76 Ge, like [31]. When the present anal- 
ysis is re-run by using the formula of [31], the T 2 ^ result is stable within 1%. 

The uncertainty related to the MaGe Monte Carlo simulation arises from two 
sources: (1) the implementation of the experimental geometry into the code (dimen- 
sions, displacements, materials); and (2) the interaction of radiation with matter (cross 
sections, final state distributions) which is performed by Geant4. These items are 
evaluated to be 1% and 2%, respectively. The estimated contribution due to the par- 
ticle tracking is based on the fact that electromagnetic physics processes provided by 
Geant4 for 7-rays and e have been systematically validated at the few-percent level 
in the energy range which is relevant for 7-ray spectroscopy [32]. In this particular 
application the Monte Carlo uncertainty is mainly due to the propagation of the ex- 
ternal 7-rays: the 2^/3/3-decay electrons generated in the germanium detectors have a 
sub-cm range and they usually deposit their entire kinetic energy, apart from small 
losses due to the escape of Bremsstrahlung or fluorescence photons. 

The combination in quadrature of the contributions reported in table 2 sums up 
to tf 2 3 %, corresponding to t°m ' 1q21 y r - 

5. Results and conclusions 

The half-life of the 2vj3f3 decay of 76 Ge was derived from the first data from the Gerda 
experiment at the INFN Gran Sasso Laboratory. Gerda operates HPGe detectors 
enriched in Ge directly immersed in liquid argon. The analysis has been carried out 
on the data collected with six enr Ge detectors (14.6 kg total weight) during 125.9 live 
days by fitting the energy spectra with a comprehensive model. The best estimate of 
the half-life of the 2i//3f3 decay is 

Tlh = (l-84^ fi t^ 6sys t)-10 21 yr = (l-84±g;ig) • 10 21 yr, (2) 
with the fit and systematic uncertainties combined in quadrature. 

The half-life is longer than all previous measurements reported in the literature. 
A summary of the Ty 2 results of 76 Ge from earlier experiments vs. publication year is 
displayed in figure 4. The figure includes nine measurements published between 1990 
and 2005 and two weighted averages. A trend towards longer Tfj 2 vames reported 
by the most recent (and lowest background) experiments is clearly seen. The Gerda 
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Figure 4. Experimental results for T^J 2 of 76 Ge vs. publication year. The 
plot includes results from the experiments ITEP-YPI [33], PNL-USC ( nat Ge) [34] 
PNL-USC-ITEP-YPI [35, 36], Heidelberg-Moscow (HdM) [11, 17] and Igex [12, 
16], as well as the re-analysis of the HdM data by Klapdor-Kleingrothaus et al 
[31] (HdM-K) and by Bakalyarov et al [37] (HdM-B). The NNDC-recommended 
value [38] and the global weighted average evaluated by Barabash [2] are also 
shown. 



result is in better agreement with the two most recent results of [31, 37] that are based 
on the re-analysis of the HdM data* . The fact that the half-lives derived in the more 
recent works — and particularly in this one — are systematically longer is probably 
related to the superior signal-to-background ratio, which lessens the relevance of the 
background modelling and subtraction. Thus the total Gerda uncertainty is compa- 
rable to what was achieved by HdM, in spite of the much smaller exposure. 
The uncertainty of the Gerda result can be further reduced in the future by accumu- 
lating more exposure and by performing a new and more precise measurement of the 
active mass of the detectors. Large data sets will reduce also the uncertainty due to 
the fit model, as background components can be better characterized and constrained 
by the (non-)observation of 7 lines. 

Using phase space factors from the improved electron wave functions reported in 
[39], the experimental matrix element for the 2vj3(3 decay of 76 Ge calculated with the 
half-life of this work is NME 2 " = 0.133±g MeV" 1 . 

The present value for NME 2l/ is 11 % smaller than that used in [5], which compares 
the matrix elements for 2v(i(i and 0^/3/3 decays. Ref. [5] also shows a relation between 
the two matrix elements for QRPA calculations of 76 Ge. According to this, the new 
value for NME 2 " results in an increase of the predicted half-life for the (V/3/3 by about 

* If T^J 2 were as short as 1.5 • 10 21 yr reported in [2], the 2f/3/3 decay of 76 Ge would account for 
nearly all counts detected in the range 600-1800 keV (expected: 8782.7, detected: 8976), thus leaving 
almost no space for any other background source in Gerda. 
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15%, which is well within the present uncertainty of the model calculation. 
The nuclear matrix elements of the 2^/3/3 decay of 76 Ge estimated from the charge 
exchange reactions (d, 2 He) and ( 3 He,t) are (0.159±0.023) MeV" 1 [7] and (0.23±0.07) 
MeV -1 [8], respectively. They both seem to be higher than the value reported in this 
work, but consistent within the uncertainties. 
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